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Abstract 

Background:: Crop plants are highly sensitive to ambient temperature, with a 1 °C difference in temperature 
sufficient to affect development and yield. Monocot crop plants are particularly vulnerable to higher temperatures 
during the reproductive and grain-filling phases. The molecular mechanisms by which temperature influences grain 
development are, however, unknown. In Arobidopsis tholiono, H2A.Z-nucleosomes coordinate transcriptional 
responses to higher temperature. We therefore investigated whether the effects of high temperature on grain 
development are mediated by H2A.Z-nucleosomes. 

Results:: We have analyzed the thermal responses of the Pooid grass, Brachypodium distachyon, a model system 
for crops. We find that H2A.Z-nucleosome occupancy is more responsive to increases in ambient temperature in 
the reproductive tissue of developing grains compared withvegetative seedlings. This difference correlates with 
strong phenotypic responses of developing grain to increased temperature, including early maturity and reduced 
yield. Conversely, temperature has limited impact on the timing of transition from the vegetative to generative 
stage, with increased temperature unable to substitute for long photoperiod induction of flowering. RNAi silencing 
of components necessary for H2A.Z-nucleosome deposition is sufficient to phenocopythe effects of warmer 
temperature on grain development. 

Conclusions:: H2A.Z-nucleosomes are important in coordinating the sensitivity of temperate grasses to increased 
temperature during grain development. Perturbing H2A.Z occupancy, through higher temperature or genetically, 
strongly reduces yield. Thus, we provide a molecular understanding of the pathways through which high 
temperature impacts on yield. These findings may be useful for breeding crops resilient to thermal stress. 
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Background 

Members of the Pooideae grass family, including wheat, 
barley, oat and rye, are a major source of human nutri- 
tion. The phenology of these crop plants, and the yield 
and quality of grain produced are significantly influ- 
enced by temperature [1,2], making them vulnerable to 
climate change [3,4]. 

The effects of temperature at various stages of cereal 
development have been extensively studied, and optimal 
temperatures determined for phenological phases from 
sowing and emergence through to grain development 
(reviewed in [5]). During vegetative stages, the effects of 
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temperature on growth are evident by the rise in leaf 
extension rates that occur as temperature increases [6,7]. 
During generative stages, the influence of temperature on 
leaf extension rate increases, suggesting that monocot 
plants have varying degrees of thermal sensitivity depend- 
ing on their developmental stage [7]. This is evident 
during late reproductive stages, where the effects of ther- 
mal stress are significantly stronger at anthesis and stages 
thereafter, compared to the double ridge stage, which is 
the earliest morphological sign of a reproductive plant [8]. 
Importantly, this includes a major effect of increasing tem- 
perature during endosperm development, with growth at 
moderately high temperatures of 27°C to 32°C reducing 
the duration of grain filling without a compensatory 
increase in the rate of grain filling, resulting in significantly 
reduced yield [9-12]. Increased temperatures also affect 
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the transcriptome of developing grain, resulting in grain at 
elevated temperatures having a more advanced develop- 
mental age [13-15]. Taken together, these results indicate 
there is a genome-wide mechanism that integrates thermal 
information into the transcriptome of developing grain. 

In Arabidopsis thaliana, H2A.Z-nucleosomes play a 
key role in mediating the effects of ambient temperature 
on the transcriptome[16]. H2A.Z-nucleosomes are fre- 
quently found at positions surrounding the transcription 
start site (TSS) [17-22]. Occupancy of H2A.Z-nucleo- 
somes at the TSS restricts access of transcriptional 
machinery into the gene body, and is reduced as tem- 
perature increases [16]. The reduced occupancy occurs 
irrespective of a given gene's transcriptional response to 
increased temperature, indicating eviction of H2A.Z is 
caused by exposure to warmer temperature and not sim- 
ply a consequence of a higher transcription rate [16]. The 
developmental phenotypes that occur when Arabidopsis 
plants are exposed to warmer temperatures, including 
accelerated flowering, are constitutively present at cooler 
temperatures in genotypes compromised in their ability 
to incorporate H2A.Z into chromatin [16,23-26]. H2A.Z- 
nucleosomes therefore provide a genome-wide mechan- 
ism by which the transcriptome can be coordinated with 
temperature to fine-tune development in response to the 
environment. 

To understand how crop plants respond to warmer tem- 
peratures we have used Br achy podium distachyon, a 
model Pooid grass and close relative of wheat and barley, 
which is a good exemplar of cereal biology and grain 
development [27,28]. We assessed the effects of tempera- 
ture on plant phenology and H2A.Z-nucleosomes of ther- 
mally responsive genes, and find that they are more 
pronounced in developing grain compared to vegetative 
seedlings. Grain from transgenic plants deficient in H2A.Z 
deposition resemble those of wild-type plants grown at 
higher temperature. Our results suggest that H2A.Z- 
nucleosomes are responsible for the increased thermal 
sensitivity of reproductive grain-filling tissue compared to 
vegetative tissue in monocot crop plants. 

Results 

Warmer ambient temperature is not sufficient to induce 
flowering in Brachypodium 

Flowering is a major phase transition in plants. In 
A. thaliana, the floral transition is highly responsive to 
ambient temperature, with growth at 27°C sufficient to 
overcome late flowering in short photoperiods [29]. Flow- 
ering in natural accessions of B, distachyon responds to a 
shift from long days (LD) to short days (SD) [30]. In our 
growth conditions, the accession Bd21 did not flower even 
after 150 days when grown in SD (14 h light/ 10 h dark) at 
22°C, indicating it is a long-day accession. Brachypodium 
also did not flower in SD at 27°C, indicating that, unlike 



Arabidopsis, increased ambient temperatures are unable to 
substitute for long day induction of flowering (Figure 1). 
The effect of increased temperature on flowering was 
further investigated in LD (20 h light/4 h dark), as well as 
after transfer from SD to LD. No acceleration of flowering 
was observed at 27°C compared to 22°C when plants were 
grown constantly in LD, with flowering at both tempera- 
tures occurring 22 days after germination (Figure la). 
Previous studies have suggested that temperate grasses are 
more thermally responsive during reproductive stages 
than vegetative development [6,7]. To test if this is the 
case in Brachypodium, we assessed the response to 
increased temperature following floral initiation. Accord- 
ingly, plants were grown at 22°C in SD until the emer- 
gence of leaf 7, and then transferred to LD for 2 days, 
which is sufficient for the induction of flowering by tran- 
scriptional activation of Flowering Locus T (Figure SI in 
Additional file 1). Thereafter, plants were either main- 
tained at 22°C or transferred to 27°C, also in LD. By limit- 
ing the plant's exposure to warmer temperatures following 
its commitment to reproductive development, head emer- 
gence was found to occur 4 days earlier at 27°C compared 
to 22°C (Figure la). These results indicate that increased 
temperatures can accelerate flowering in the presence of 
favorable photoperiods, but are not alone sufficient to trig- 
ger flowering. 

This lack of responsiveness of vegetative plants to 27°C 
could be because Brachypodium has evolved to respond to 
a different range of temperatures than Arabidopsis, or 
because it is differentially responsive to temperature 
depending on its developmental phase. To assess the latter 
possibility, we assayed the effects of high temperature on 
grain development, because it is a late-reproductive stage 
that is particularly sensitive to thermal stress [2,11]. Plants 
were grown at 22/17°C day/night cycles until the onset of 
endosperm development [28], and then maintained at this 
temperature or shifted to 27/22°C. These temperature 
treatments were used because growth between 27 and 
30°C is recognized as a moderately high thermal stress 
that adversely affects yield in wheat and barley, compared 
to control conditions at 20 to 24°C [11,12,15,31]. At 27°C, 
grain weight initially increased faster than at 22°C, such 
that after 4 days at 27°C, grain was 14% heavier than at 
22°C {P < 0.005), indicating an accelerated rate of grain 
development (Figure 2a,b). However, grain-filling ceased 
10 days after pollination (DAP) at 27°C while it continued 
until 16 DAP at 22°C, which caused a significant reduction 
in final fresh and dry weights for grain grown at 27°C. The 
reduced duration of grain-filling was confirmed with 
mature grain weight measurements that show a 16% 
decrease {P < 0.001) at the higher temperature (Figure 2c). 
Therefore, grain development is significantly affected by 
increased ambient temperature, with higher temperatures 
accelerating the rate of grain development so that there is 
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Figure 1 Higher ambient temperature is not sufficient to induce flowering in Brachypodium distachyon . (a) Plants were grown in one of 
three photoperiod conditions: SD (14 h light/10 h dark), LD (20 h/4 h) or in LD after shifting from SD. In each condition, plants were grown at 
either 22°C (grey) or 27°C (black). DNF (did not flower) indicates the non-flowering phenotype of plants grown in SD for 150 days. Values are the 
mean ± standard error often plants. {**''P < 0.001). (b) SD grown plants after 50 days at constant temperature of 22°C and 27°C. 



a reduced developmental window during which grain fill 
occurs, causing a reduction in yield. 

The Sracfjypod/umtranscriptome responds to changes in 
ambient temperature 

Having observed that reproductive tissue appears to be 
more developmentally responsive to increased tempera- 
ture than vegetative plants, we sought to test if this is 
reflected in the regulation of transcription and chroma- 
tin state by temperature. To identify marker genes up- 
regulated in response to increased ambient temperature, 
we assayed the transcriptional responses of vegetative 
plants shifted from 12°C to 22°C or 27°C for 2 and 
24 hours, using whole genome microarrays and quantita- 
tive real-time PGR (qRT-PCR), to obtain rapid and longer- 
term transcriptomic responses to temperature change. In 
Brachypodium seedlings, 9% of all expressed transcripts 
were either significantly up- or down-regulated at least 



two-fold after the shift from 12°C in at least one of the 
temperature treatments (two-way ANOVA, temperature 
effect P < 0.05; Figure 3a-c; Figure S2a-b in Additional file 
1). Importantly, as in other systems [32], transfer to 27°C 
is below the threshold for inducing most heat shock genes, 
while HSP70, HSF23 and HSP90 are strongly induced 
(Figure 3d). Bradi4g32941 and Bradilg32990 were also 
selected as genes for further analysis as they too were up- 
regulated at 27°C (Figure 3e). This induction occurs both 
in response to short- and medium-term shifts as well as in 
plants grown continuously at the warmer temperature, 
suggesting that these genes are responding to absolute 
temperature, and are not dependent on a temperature 
change for their expression to be induced (Figure 3e,h). 
We also identified genes whose transcript levels were 
down-regulated, such as Bradi2gl4220, BradiSg00970, 
Bradi2g484S0 and Bradi4gl7230 (Figure 3f), as well as 
genes that were unaffected by increased temperature. 
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Figure 2 Increased ambient temperature reduces grain yield in Brachypodium. {a,b) Fresh weight (a) and dry weight (b) of developing 
grain from plants grown at either 22/1 7°C (black line) or 27/22°C (red line) (day/night temperatures). Data are the replicate of 5 replicate plants, 
including measurements of at least 20 seed. Values are the mean ± standard error (**P < 0.01). (c) Final yield measurements in units of weight 
per 10 grain for plants transferred from 22/1 7°C to 27/22°C throughout grain-filling, then transferred back to 22/1 7°C at 16 days after pollination, 
compared to plants grown constantly at 22/1 7°C (***P < 0.001). Data are from 20 biological repeats. Values are the mean ± standard error. 
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Figure 3 The 5rac/7ypoc//umtranscriptome responds to changes in ambient temperature, (a) Transcript profiling experiment shows a 
robust response to changes in ambient temperature in vegetative seedlings. The heat map depicts all differentially expressed genes (DEGs) with 
at least two-fold change in any of the temperature treatments as determined from two-way ANOVA {P for temperature effect < 0.05). Expression 
levels of up-regulated genes are in shades of red and of down-regulated genes in shades of green. (b,c) Venn diagrams of total numbers of up- 
regulated (b) or down-regulated (c) DEGs in vegetative seedlings after 24 h shift to either 22°C (green), 27°C (red), or in both temperatures 
(yellow). The two-tailed P-values for the significance of the overlap represented under the Venn diagrams have been calculated using Fisher's 
exact test. {6)HSF23 (green line), H5P70 (black line) and HSP90 (red line) are induced strongly with increasing temperature, in contrast to other 
heat-shock genes (grey) that do not respond significantly over the temperature range assessed, (e-g) Quantitative real-time PGR (qRT-PGR) 
analysis of genes that are up-regulated by increasing temperature (e), down-regulated (f) or show constant expression (g) in vegetative seedlings 
24 h after temperature-shift. (h)qRT-PGR analysis of up-regulated genes in plants grown constantly at either 17°G, 22°G or 27°G. (iJ)qRT-PGR 
analysis of genes in developing grain that are up-regulated by increasing temperature (i), or remain constant within the temperature range (j). 
(k)qRT-PGR analysis of genes that are up-regulated by temperature with known roles in developing grain. Data are from at least three biological 
replicates. 
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including Bradi3g31120 and Bradilg47790 (Figure 3g). 
These results indicate that while the higher temperature 
regime does not trigger flowering in short photoperiods, 
the transcriptome of young seedlings remains thermally 
responsive. 

To determine if this transcriptional response to tem- 
perature is tissue specific, we analyzed the expression of 
these marker genes in developing grain grown at either 
22/17°C (day/night) or 27/22°C. We assessed the tran- 
scriptional response during endosperm development as 
this is the stage of grain development in wheat that 
responds strongly to increased temperature [11,15]. 
Hence, at 6 to 8 DAP, plants were either shifted to 27/ 
22°C or maintained at 22/17°C, and grains were har- 
vested after 24 h. We found that HSF23, HSP90, Bra- 
di4g32941 and Bradilg32990 were also up-regulated in 
developing grain (Figure 3i) and that Bradi3g31120 and 
Bradilg47790 transcript levels remained constant 
between the two temperature treatments (Figure 3j). 
Expression levels of HSP70 were very low and did not 
appear to be thermally regulated in developing grain, 
which is similar to results in Arabidopsis where HSP70 
was not temperature-regulated in seeds [33]. We could 
not detect expression of the genes that weredown-regu- 
lated by higher temperature in young seedlings. We also 
assessed the expression of some genes with important 
roles during grain development (Figure 3k) shown to be 
thermally responsive during grain fiUing in wheat [13], 
including beta-amylase {AMYl) and UDP-glucose pyro- 
phosphorylase {UDP-GPP) from the starch metabolism 
pathway, and Serpin 2A, which has a role in plant 
defense [13]. These genes were up-regulated in grain 
grown at 27/22°C, relative to 22/17°C, which is consistent 
with observations in wheat. These results indicate that 
the developing grain of Brachypodium responds tran- 
scriptionally to increased ambient temperature, and that 
some genes display a shared transcriptional response to 
thermal stress in young seedlings and developing grain. 

H2A.Z-nucleosomes are more thermally responsive in 
developing grain than vegetative seedlings 

The histone variant H2A.Z has an important role in regu- 
lating the ambient temperature transcriptome of Arabidop- 
sis, and loss of H2A.Z incorporation into nucleosomes 
surrounding TSSs promotes thermally sensitive develop- 
mental responses [16,26]. We therefore investigated the 
behavior of H2A.Z-nucleosomes in young seedlings and 
developing grain to determine if these nucleosomes might 
explain the altered developmental responsiveness of the 
two tissue types to increased temperature. Nucleosome 
positions surrounding the TSS were predicted using soft- 
ware prediction algorithms and confirmed using micrococ- 
cal nuclease (MNase) digestion and subsequent histone 3 
(H3) chromatin-immunoprecipitation (ChIP) (Figure 4a,b; 



Figure S3 in Additional file 1). For mF23 and HSP70, we 
used quantitative-PCR (qPCR) with primer pairs tiled 
across the promoter region surrounding the TSS, which 
revealed the appearance of two characteristic peaks that 
are suggestive of -1 and +1 nucleosomes (Figure 4a,b). As 
the software prediction algorithms accurately predicted 
nucleosome positions for HSF23 and HSP70, we used the 
software to identify the likely +1 nucleosomes of the other 
thermally responsive genes, which were confirmed by 
qPCR on DNA isolated from H3 ChIP experiments [34] 
(Figure 4c-e). 

To investigate if H2A.Z is incorporated into these 
nucleosomes, we surveyed the Brachypodium genome for 
H2A.Z homologues. Phylogenetic analysis revealed the pre- 
sence of three H2A.Z genes in Brachypodium: BdHTAl, 
BdHTA9 and BdHTAll (Figure 5a; Figure S4 in Additional 
file 1). Transcript analysis in different tissue types, includ- 
ing young seedlings and developing grain, showed that 
BdHTAl and BdHTA9 are strongly expressed, while we 
could not detect transcripts of BdHTAll (Figure 5b). We 
therefore selected BdHTA9, based on its homology to the 
H2A.Z orthologue investigated in Arabidopsis {AtHTAll) 
and its expression in seedlings and developing grain, for 
further analysis. To study the dynamics of BdHTA9 in 
response to temperature, we created a BdHTA9 3XFLAG 
form expressed under its native promoter. ChIP analysis 
using anti-FLAG resin on Brachypodium seedlings grown 
at 17°C revealed the presence of HTA9 at the presumptive 
-1 and +1 nucleosomes of both HSF23 and HSP70, and +1 
nucleosomes of the other thermally responsive genes 
(Figure 5c,d). These results are consistent with localization 
studies of H2A.Z in other organisms, which show incor- 
poration into nucleosomes surrounding the TSS [17-22]. 

To determine if eviction of H2A.Z-nucleosomes occurs 
in vegetative seedlings, as has been described in Arabidop- 
5/5 [16], we performed nucleosome analysis and ChIP for 
BdHTA9 3X FLAG from plants maintained at 17°C or 
shifted to 27°C. For all genes and nucleosome sites ana- 
lyzed, we did not detect a decrease in H2A.Z and nucleo- 
some occupancy between 17°C and 27°C (Figure 6a-c; 
Figure S5 in Additional file 1). This appears to be indepen- 
dent of transcriptional status, since it was observed equally 
for genes up-regulated, down-regulated or with constant 
expression after shifting to 27°C (Figure 6a-c). 

This apparent stability of H2A.Z-nucleosomes to higher 
temperature in vegetative plants could be a general feature 
of chromatin in Brachypodium, or it could be an indicator 
of temperature responsiveness being dependent on devel- 
opmental phase. To test this latter possibility, we exam- 
ined the thermal responsiveness of H2A.Z-nucleosomes in 
developing grain. To ensure that any altered molecular 
response to temperature was not a consequence of a dif- 
ference in the physical temperature of the two tissue types, 
thermal imaging was obtained for both organs in each 
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Figure 4 Nucleosome positioning by H3 ChIP analysis. (a,b)ChlP of cross-linked H3 at the promoter sites of HSF23 (a) and HSP70 (b) reveal 
well-positioned -1 and +1 nucleosomes. The x-axis indicates the central position of each amplicon relative to the TSS. In each schematic, the 
promoter (solid line), 5'UTR (white box), exons (black box) and TSS (arrow) are shown. (c-e)ChlP of cross-linked H3 for genes that were detected 
to be up-regulated (c), constant (d) or down-regulated (e) in response to temperature reveal sites that show strong enrichment of H3 at 
probable +1 nucleosome sites according to in s/7/cosequence analysis (see Materials and methods). The x-axis indicates the central position of 
each amplicon relative to the TSS. Values from H3 and mock reactions are shown in black and grey, respectively. Values are the mean ± 
standard error of three biological replicates. 



treatment (Figure S6 in Additional file 1). This analysis 
showed that there was no difference in physical tempera- 
ture between vegetative seedlings and developing grain 
(Figure S6 in Additional file 1). We then investigated H2A. 
Z-nucleosome behavior under the same conditions used 
for transcript analysis. BdHTA9 3X FLAG plants were 
grown at 22/17°C until 6 DAP, and then shifted to 27/22°C 
or maintained at 22/17°C, with grains harvested after 24 h. 
Chip and nucleosome analysis showed that there is a strik- 
ing reduction in HTA9 occupancy at 27/22°C, relative to 
22/17°C, for all genes analyzed, including AMW, UDP-GPP 
and Serpin 2A (Figure 6d-f). We also observed a compar- 
able reduction in nucleosome levels (Figure S7 in Addi- 
tional file 1). These results show that in the endosperm of 
developing grain there is a considerable increase in the 
mobility of H2A.Z at the +1 nucleosomes of these genes. 



compared to the same nucleosome sites in vegetative seed- 
lings. In developing grain the transcript levels of BdHTA9 
did not change with temperature, indicating that the 
increased mobility of H2A.Z in grain does not occur as a 
consequence of reduced BdHTA9 transcription (Figure S8 
in Additional file 1). This result is consistent with the 
absence of large-scale phenological changes in vegetative 
plants at 27°C, while a considerable acceleration in grain 
development is seen at this higher temperature. 

Grain from ARP6 RNAi lines phenocopy responses to 
increased temperature 

To determine whether loss of H2A.Z-nucleosomes at higher 
temperature is the cause of altered grain development, or 
simply a correlated event, we sought to perturb H2A.Z- 
nucleosome occupancy independently of temperature. H2A. 
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Figure 5 Identification and nucleosome positioning of 8cfHTA9 in Brachy podium, (a) An unrooted maximum-likelihood phylogenetic tree 
of HTA proteins in Brochypodium, Arobidopsis, humans and yeast, constructed using MEGA5 with 100 bootstrap replicates, summarizes the 
evolutionary relationship among the HTA proteins and the separation in four phylogenetic subfamilies. Branches are drawn to scale with scale 
bar representing the number of substitutions per site, (b) Reverse-transcriptase PGR analysis of BdHTAI, BdHTA9 and BdHTAll from leaf, apex and 
endosperm tissue of Bd21. Two biological replicates are shown. (c,d)ChlP analysis of HTA9:3XFLAG (H2A.Z) at 17°C shows H2A.Z is enriched at 
the -1 and +1 nucleosomes of HSF23 (c) and HSP70 (d). The x-axis indicates the central position of each amplicon relative to the TSS. In each 
schematic, the promoter (solid line), 5'UTR (white box), exons (black box) and TSS (arrow) are shown. Mock reactions (grey) were performed on 
identical tissue from wild-type plants. Values are the mean ± standard error of three biological replicates. 



Z is inserted into chromatin through the highly conserved 
SWRl complex [24,35,36]. We disrupted the activity of 
SWRl using RNA interference (RNAi) silencing of a key 
conserved component, ACTIN RELATED PROTEIN 6 
{ARP6), of which there is only one copy in the Br achy po- 
dium genome {Bradi2gl0130), Using three independent 
transgenic lines with reduced expression of ARP6, we 
observed a reduction in seed weight (60%, P < 0.001) and 
overall yield per plant (30 to 50%, P < 0.01) that was com- 
parable with the reduced yield observed in wild-type plants 
grown at 27/22°C (Figure 7b, c). We also observed an 
increased rate of floret sterility and grain abortion in the 
ARP6 RNAi lines (Figure 7c; Table SI in Additional file 1), 
which is consistent with observations in wheat where high 



temperature treatment of early developing grain, prior to the 
onset of endosperm development, promotes grain abortion 
and reduction in grain number [37-39]. Other than the phe- 
notypes associated with grain development, the ARP6 RNAi 
lines appeared phenotypically normal and displayed a flow- 
ering-time consistent with that of wild-type plants, including 
an inability to flower in SD (Figure S9 in Additional file 1). 
These results show that developing grain oiARP6 RNAi 
plants grown at 22/17°C simulate grain harvested fi:om wild- 
type plants grown at 27/22°C. 

To determine if the grain development phenotypes 
observed in ARP6 RNAi plants may be caused by the 
mis-regulated expression of thermally responsive genes, 
we compared the transcript levels of thermally induced 
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Figure 6 Occupancy of H2A.Z-nucleosomes is reduced at higher ambient temperatures in developing grain but not in vegetative 
seedlings. (a-c)ChlP analysis of HTA9:3XFLAG (H2A.Z) at 22°C and 27°C in vegetative seedlings at +1 nucleosomes of genes whose expression 
was up-regulated (a), remained constant (b) or down-regulated (c) upon an increase in temperature. (d-f)ChlP analysis of HTA9:3XFLAG (H2A.Z) 
at 22/1 7°C and 22/27°C in developing grain at +1 nucleosomes of genes whose expression was up-regulated (d) in both seedlings and grain, or 
remained constant (e) upon an increase in temperature. (f)ChlP analysis of HTA9:3XFLAG (H2A.Z) at 22/1 7°C and 27/22°C for genes with roles in 
grain development. Mock reactions (grey) were performed on identical tissue from wild-type plants. < 0.001. 



genes in wild-type and transgenic plants grown at 22/17°C. 
We used qRT-PCR to assess the relative expression of 
genes that were up-regulated in developing grain upon 
transition from 22/17°C to 27/22°C. These results show 
that reduced activity of ARP6 causes elevated transcript 
levels for genes that were induced by increased tempera- 
ture in wild-type plants (Figure 8). Some variation in the 
behavior of transcript levels between RNAi lines may 
reflect residual ARP6 activity. These results indicate that 
when measured transcriptionally, the developing grain of 
ARP6 RNAi lines grown at 22/17°C resemble those of 
wild-type plants grown at warmer temperatures. 

Discussion 

Ambient temperature has a large effect on plant devel- 
opment, reflected in changes to the transcriptome 
[32,40]. In A. thaliana, H2A.Z-containing nucleosomes 



are necessary to correctly express the temperature tran- 
scriptome[16]. We therefore sought to determine 
whether H2A.Z-mediated perception of ambient tem- 
perature is conserved in monocot plants, and whether 
the response to increased temperature varies between 
developing grain and vegetative seedlings. 

We have observed that chromatin sites containing 
H2A.Z-nucleosomes are more responsive to increases in 
ambient temperature in developing grains than vegeta- 
tive seedlings. These molecular observations support 
phenological changes reported in this and other studies, 
which show that increased temperatures cause a greater 
response in plants progressing through late-reproductive 
development compared to those of vegetative stages. For 
example, the base temperature used to calculate the rate 
of development using thermal time in wheat is not con- 
stant throughout plant development, with the base 
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Figure 7 Seeds from plants with reduced expression of ARP6 phenocopy seed from plants grown at higher temperature, (a) Relative 
expression of ARP6 in wild-type {WT; Bd21) and three independent ARP6 RNAi transgenic lines (Ji generation); < 0.001. (b,c) Seed weight 
measurements (b) and yield per plant (c) in WT (Bd21) at 22/1 7°C and 27/22°C, as well as three independent transgenic lines at 22/1 7°C (Ti 
generation). Data are the mean ± standard error of at least 15 grains < 0.01; **''P < 0.001). (d) A representative spike from WT Bd21 and 
ARP6 RNAi.1 displaying the empty florets (white arrowheads) that contained aborted grain. Scale bar, 1 cm. 
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Figure 8 Genes up-regulated by increased ambient temperature are up-regulated in seeds of ARP6 /?AM/ transgenic lines at 22°C. Expression 
of genes that are up-regulated in developing grain by transfer from 22/1 7°C to 27/22°C (Figure 3) are up-regulated in grain of three ARP6 RNAi 
transgenic lines (Ti generation) grown at 22/1 7°C. Values are the mean ± standard error of 4 replicates, each containing 6 grain at 8 DAP. 
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temperatures for terminal spikelet to anthesis, and 
anthesis to maturity being considerably higher than that 
of the vegetative stages [41-43]. A differential response 
of vegetative and reproductive tissue to moderate and 
high temperature has also been observed in other plants, 
such as orchids and strawberry [44,45]. Similarly, it has 
been shown in Arabidopsis that cold temperature treat- 
ment causes approximately ten-fold more genes to be 
induced in seeds than in seedlings [33]. One hypothesis 
for the distinct responses of vegetative and reproductive 
tissues is that it reflects an adaptation to low tempera- 
tures during early developmental stages and to warmer 
conditions during reproductive stages [42,43]. Our findings 
suggest that plants may be able to modulate responsiveness 
to temperature in a tissue- and lifecycle-dependent manner 
by controlling the stability of H2A.Z-nucleosomes. To 
determine how direct the link between H2A.Z-nucleo- 
somes is to the control of grain development in response 
to temperature, it will be necessary to identify the loci reg- 
ulating this process and the binding dynamics of H2A.Z- 
nucleosomes to their promoters. 

The effect of temperature on floral initiation has been an 
active area of research in plant development. Here we find 
that increased temperature alone cannot promote flower- 
ing in non-inductive photoperiods, indicating that Brachy- 
podium is an obligate long-day plant and that increased 
ambient temperature cannot compensate for photoperiod 
as a floral inductive signal. These observations are consis- 
tent with other studies, where it has been shown in wheat 
and barley that increased temperature in SD conditions 
does not promote floral initiation, but in fact delays the 
onset of reproductive development [46,47]. Similarly, in 
strawberry, olive, rice flowers, perennial grasses and orch- 
ids moderately high temperatures delay or completely 
inhibit floral induction [44,45,48-50]. Taken together, 
these observations suggest that some plants, including 
annual temperate grasses such as Br achy podium, have 
adopted different strategies to respond to increased ambi- 
ent temperature compared to Arabidopsis, where flower- 
ing is accelerated by high temperature. 

The results presented in this study indicate that the 
transcriptional responses of ripening grain to increased 
temperature, and the accompanied accelerated rate of 
development and reduction in grain weight, are at least 
partially mediated by H2A.Z-nucleosomes. Recent studies 
in wheat and barley have investigated the transcriptional 
events initiated by moderate and high increases in tem- 
perature during grain development [13-15]. These studies 
have shown that the transcriptome of a developing grain is 
highly responsive to increased ambient temperatures [14] 
and that the developmental age of the grain is more 
advanced at warmer temperatures, as measured transcrip- 
tionally [15]. Therefore, it may be possible to reduce the 
negative effects of thermal stress on grain weight and 



quality by nullifying the transcriptional response of devel- 
oping grain to increasing temperatures. This is a concept 
that has already been investigated in rice, where attenuat- 
ing the transcriptional activity of a -amylases, which are 
normally up-regulated by increased temperatures, has 
been used to ameliorate the deleterious effects that high 
temperatures have on grain quality [51]. Given the broad 
role for H2A.Z-nucleosomes in regulating the ambient 
temperature transcriptome in Arabidopsis seedlings, it is 
possible that a large proportion of the transcriptional 
responses initiated by increased temperature in developing 
grain are coordinated by H2A.Z-nucleosomes. Therefore, 
reducing the responsiveness of H2A.Z-containing nucleo- 
somes in developing grain may help improve yield and 
grain quality for plants grown in high temperatures. This 
hypothesis is supported by findings in Arabidopsis that 
plants which lack H2A.Z have fertility defects [26], as well 
as the results of this study that decreased activity of ARP6 
causes reduced seed weight and increased expression of 
thermally responsive genes in plants grown at lower 
temperatures. 

A key question for future research will be to understand 
the mechanism underlying the differential response of 
H2A.Z-nucleosomes to increased temperature in reproduc- 
tive and vegetative tissue. Histone modification is one pos- 
sible explanation for the more pronounced effect of 
thermal stress on H2A.Z behavior in nucleosomes of devel- 
oping grain compared to vegetative seedlings. Recent 
results in humans have shown that while H2A.Z is loca- 
lized at sites surrounding the TSS of both active and poised 
gene promoters, acetylated H2A.Z is only found at TSSs of 
active genes [52]. In yeast, acetylation of H2A.Z was found 
to correlate with genome- wide gene activity [53], and con- 
fer with nucleosome destabilization and an open conforma- 
tion of chromatin [54]. It is possible, therefore, that 
differences in the degree of H2A.Z acetylation may account 
for tissue-specific variation in thermal responsiveness. In 
addition to the H2A.Z-nucleosome dependent regulation 
of transcription, we also observe warm temperature induc- 
tion of gene expression in seedlings, where no significant 
decreases in H2A.Z-nucleosomes are seen. This is consis- 
tent with work in Arabidopsis showing that H2A.Z-nucleo- 
somes do not account for all the transcriptional responses 
to warmer temperature [16]. The identification of these 
parallel temperature responsive pathways will be of 
interest. 

The results presented in this work, as well as many phe- 
nological studies, suggest that Brachypodium and other 
grasses are more responsive to warmer temperatures as 
they enter the generative phase of their lifecycle. This may 
be because the initiation of flowering represents a point of 
no return, since the reproductive program must be com- 
pleted. Higher temperatures would then represent an 
environmental cue associated with the likely onset of 
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temperature stress and drought. Accelerating the develop- 
mental transition through this vulnerable stage to the resili- 
ent desiccated grain state may therefore be an adaptive 
response to anticipated abiotic stress. Conversely, activating 
flowering in response to warmer temperatures indepen- 
dently of photoperiod would risk exposing the inflores- 
cence and the developing grains to frost [55]. This study 
highlights the value of complementing work in Arabidopsis 
with monocot systems to assess the direct applicability of 
Arabidopsis research to other systems, especially for traits 
that are important for crop yield, such as temperature 
perception. 

Conclusions 

In common with other monocot plants, Brachypodium 
displays a marked difference in temperature responsive- 
ness depending on the developmental stage. We show that 
the process of grain-fill is particularly sensitive to warmer 
temperatures. Our data indicate that H2A.Z-nucleosomes 
are necessary for coordinating the sensitivity of Brachypo- 
dium to increased temperature during grain development. 
Perturbing H2A.Z occupancy, through higher temperature 
or genetically, strongly reduces yield. Thus, we provide a 
molecular understanding of the pathways through which 
high temperature impacts on yield. These findings may be 
useful for breeding crops resilient to thermal stress. 

Materials and methods 

Plant material and growth conditions 

All experiments were carried out using Brachypodium dis- 
tachyon . Material included wild-type Bd21 and trans- 
formed lines containing BdHTA9 with a carboxy-terminal 
3X FLAG tag, including native promoter with all exons 
and introns, and transformed lines containing a RNAi 
construct containing sequence complementary to ACTIN 
RELATED PROTEIN 6 {ARP6) (nucleotides 360 to 508 of 
open reading frame). Transformation was performed as 
described in [56]. All plants were grown in growth cabi- 
nets (Microclima MClOOOE, Snijders Scientific, Tilburg, 
Netherlands) with 70% relative air humidity and 350 (imol 
m'^ s'^ PPFD (Photosynthetic Photon Flux Density). 

For flowering-time experiments, plants were grown in 
one of three photoperiod regimes: i) SD (14 h light/10 h 
dark); ii) LD (20 h/4 h); or iii) transferred from SD to LD 
at the seven leaf stage. Plants were grown at constant 
temperatures of either 22°C or 27°C. For analysis of tem- 
perature effects on grain development, plants were grown 
in SD at 22/17°C (day/night temperatures) until the 
seven-leaf stage to synchronize development, and shifted 
to LD at the same temperature. After head emergence, 
date of pollination was monitored for the main stem. At 
5 DAP, plants were either maintained at 22/17°C or 
shifted to 27/22°C and developing grains were collected 



at two-day intervals. At 16 DAP, plants grown at 27/22°C 
were shifted back to 22/17°C until maturity to measure 
final yield. Dry weight of developing grain was measured 
after desiccating grain at 60°C for 10 days. For each time 
point and temperature, a total of 20 seed were measured 
that were collected from 5 plants. 

For transcript analysis of vegetative plants, seedlings 
were grown at 17°C in SD until emergence of the third 
leaf, and transferred to 12°C for 48 h before shifting to 
either 22°C or 27°C. Samples were collected at 2 h and 
24 h post-shift. Transcript analysis of developing grain was 
performed only on grain harvested from the main stem. 
Grain were collected from plants grown at 22/17°C until 
the end of 6 DAP, which were then transferred to 27/22°C 
or maintained at 22/17°C. Samples were collected 24 h 
after shift. Each replicate contained pools of 6 grain, from 
which the glumes had been removed. 

Transcript analysis 

Transcript analysis was performed on total RNA extracted 
using the Plant RNeasy Plant Mini kit (Qiagen, Crawley, 
West Sussex, UK). Single-stranded cDNA was amplified 
for microarray analysis using the Ambion WT Expression 
Kit (Ambion Life Technologies, Paisley, Renfrewshire, 
UK), labeled and hybridized to the custom Brachypodium 
microarray (Affymetrix, Santa Clara, California, USA) 
using the GeneChip WT terminal labeling and hybridiza- 
tion kit (Affymetrix). Microarray data were analyzed using 
GeneSpringGX v.ll (Agilent, Santa Clara, California, 
USA). For each time-point after the temperature shift (2 h 
and 24 h) and temperature treatment (12°C, 22°C and 
27°C), two to three biological replicates were analyzed. 
Samples taken before the shift at 12°C 0 h were considered 
to be equivalent to the 12°C 24 h treatment for the pur- 
pose of the statistical analysis. Raw data were normalized 
using the RMA algorithm and filtered by expression level 
using the cut-off value of >20 in at least one of the treat- 
ments to define the group of genes expressed in vegetative 
seedlings. The statistical significance of changes in the 
transcripts with expression levels above the cut-off level 
was assessed using two-way ANOVA with temperature 
and time as main factors. The P-values were adjusted for 
multiple testing by the Benjamini and Hochberg s method 
at a false discovery rate of <0.05. The transcripts signifi- 
cantly affected by temperature as the main factor (cor- 
rected P < 0.05) and displaying >2-fold change in any of 
the temperature-time combinations were reported as dif- 
ferentially expressed genes (DEGs). Hierarchical clustering 
and heat mapping were used to visually display DEGs. The 
significance of the overlap between the lists of DEGs in 
different temperature-time treatments represented in 
Venn diagrams were tested using Fisher's exact test. The 
DEGs were annotated using the bradinet platform of web 
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tools [57], and are provided as a list in Table S5 in Addi- 
tional file 2). The data for these experiments have been 
deposited in ArrayExpress (E-MEXP-3918). 

For transcript analysis by qRT-PCR, cDNA was synthe- 
sized using the Fermentas First Strand cDNA synthesis kit 
(Thermo Scientific, St. Leon-Rot, Germany). qRT-PCR 
was performed using SYBR Green I (Roche, Burgess Hill, 
West Sussex, UK) and a Roche LightCycler 480. All qRT- 
PCR data points are the average of three biological repli- 
cates, with two technical replicates performed in each 
reaction. See Table S2 in Additional file 3 for oligonucleo- 
tide sequences used for qRT-PCR. Expression of candidate 
genes was normalized against SamDC and UBC18[58], 

Nucleosome positioning and ChIP 

Nucleosome positioning and ChIP were performed 
essentially as described [16], with minor modifications. 
For positioning of nucleosomes in vegetative plants, 
chromatin from seedlings grown at 17°C was cross-linked 
using 1% formaldehyde. Chromatin was fragmented using 
0.2 units of micrococcal nuclease (Sigma, Gillingham, 
Dorset, UK) as described [16]. Nucleosome positions 
were identified by immunopurification using the H3 anti- 
body (Abeam, abl791, Cambridge, Cambridgeshire, 
UK)). Relative enrichment of associated DNA fragments 
was analyzed by qPCR. Nucleosome occupancy was 
determined as a fraction of uncut chromatin DNA, with 
data plotted against the HSF23 and HSP70 promoters. 
HSP70 was selected for detailed nucleosome positioning 
because it is the homologous gene used for transcript 
and H2A.Z-nucleosome analysis in A. thaliana[16], and 
HSF23 was selected because it is expressed at a level pro- 
portionate to ambient temperature within the range 12 
to 27°C. For in silico predictions of +1 nucleosome sites, 
1 kb and 250 bp of sequence upstream and downstream 
of the TSS, respectively, were used to query the online 
nucleosome position prediction software (version 3.0) on 
the Segal lab website [34]. The two most probable sites 
were selected for qPCR analysis of H3 ChIP DNA, with 
amplicons that provided the strongest signal used for 
further analysis. The center of the amplicon indicates its 
position relative to the TSS. Oligonucleotide sequences 
are provided in Table S3 in Additional file 4. 

Similar methods as outlined above were used to deter- 
mine nucleosome and H2A.Z dynamics at different tem- 
peratures. Nucleosome and H2A.Z dynamics were 
assessed using chromatin isolated from seedlings that were 
grown at 17°C until 10 days after germination, and then 
either maintained at 17°C or shifted to 27°C for 24 h. 
Chromatin was digested with MNase, and for nucleosome 
analysis, mononucleosome-sized fragments were gel puri- 
fied and used in qPCR. For H2A.Z dynamics, ChIP was 
performed on MNase digested chromatin using the 
ANTI-FLAG M2 Affinity Gel (Sigma), and histone-DNA 



complexes eluted using the 3X Flag peptide (Sigma). For 
assays of nucleosome and H2A.Z occupancy in developing 
grain, plants were grown at 22/17°C until 6 DAP. Plants 
were then either maintained at 22/17''C or shifted to 27/ 
22°C and developing grains (1 g per ChIP assay) were col- 
lected after 24 h (7 DAP), with grain harvested only from 
the main stem. Native ChIP assays were performed as 
described above, with the following modifications. Ground 
tissue was washed twice in extraction buffer I, and three 
times in extraction buffer II to remove metabolic contami- 
nants. Following MNase digestion, two rounds of centrifu- 
gation were performed on uncut chromatin DNA and 
digested chromatin DNA prior to reverse cross-linking 
and immunopurification, respectively, to further remove 
contaminants. Nucleosome and H2A.Z occupancy was 
determined as a fraction of input DNA that was diluted 
ten-fold for qRT-PCR. Oligonucleotide sequences are pro- 
vided in Table S3 in Additional file 4. All ChIP assays 
were performed at least three times, and two technical 
replicates were performed for each qRT-PCR. 

Sequence alignment and phylogenetic analysis 

Brachypodium HTA sequences were identified by 
BLAST search [59] from public databases using Arabi- 
dopsis HTA protein sequences as bait (Table S4 in 
Additional file 5). Multiple sequence alignment was per- 
formed using webPRANK[60], and manually corrected 
using GENEDOC [61]. Phylogenetic trees were con- 
structed in MEGA5 [62] using the maximum likelihood 
algorithm with the following parameters: JTT substitu- 
tion model. Gamma distributed rates and partial dele- 
tions. Bootstrap values are based on 100 replicates for 
testing the significance of the nodes. 

Seed weight measurements from ARP6 RNAi lines 

For wild-type (Bd21) and each transgenic ARP6 RNAi 
line, seed weight measurements were determined for 15 
seeds harvested from 3 plants (5 seeds from each repli- 
cate plant). Seeds were harvested from basal florets of 
spikelets from the main spike, and the lemma and palea 
were removed before weighing. For yield per plant mea- 
surements, all seeds were harvested from three replicate 
plants. 

Statistical analysis 

When not described specifically, differences between 
treatments were tested by Student's ^-test. Results in fig- 
ures are shown as means ± standard error. 

Additional material 



Additional file 1: Supplementary figures. 

Additional data file 2: Table 55. The differentially expressed ger^es 
(DEGs), defir^ed as significantly affected by temperature as the main 
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factor {P < 0.05) and displaying >2-fold change in at least one of the 
temperature-time treatments. 

Additional file 3: Table S2. 

Additional file 4: Table S3. 

Additional file 5: Table S4. 

Additional file 6: Supplementary methods and figure legends. 
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